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ABSTRACT 

Recent observations of the black hole (BH) - bulge scaling relations usually report positive redshift evolution, 
with higher redshift galaxies harboring more massive BHs than expected from the local relations. All of these 
studies focus on broad line quasars with BH mass estimated from virial estimators based on single-epoch 
spectra. Since the sample selection is largely based on quasar luminosity, the cosmic scatter in the BH-bulge 
relation introduces a statis tical bias leading to on average more massive BHs given galaxy properties at high 
redshift dLauer et alj 2007). We here emphasize a previously under-appreciated statistical bias resulting from 
the uncertainty of single-epoch virial BH mass estimators and the shape of the underlying (true ) BH mass 
functio n, which leads to on average overestimation of the true BH masses at the high-mass end (IShen et alJ 
2008b). We demonstrate that the latter virial mass bias can contribute a substantial amount to the observed 
excess in BH mass at fixed bulge properties, comparable to the Lauer et al. bias. The virial mass bias is 
independent of the Lauer et al. bias, hence if both biases are at work, they can largely (or even fully) account 
for the observed BH mass excess at high redshift. 

Subject headings: black hole physics — galaxies: active — quasars: general — surveys 



1. INTRODUCTION 

The redshift evolution of the local scaling relations between 
galaxy bulge properties and the mass of the central super- 
massive black hole (SMBH) has important clues to the es- 
tablishment of these tight relations across cosmic time. There 
is currently a huge effort in measuring this evolution, either 
in terms of the BH mass-bulge velocity dispersion (M. - 
er) relation, or the BH mass-bulge stel lar mass/luminosity 
(M, - Mbuipe. M, - Lbuiee) relation (e.g., IShields et al] 12003c 

""20081 



Peng et ail 



Woo et all 



2006ab; Salviander et al. 2007; Shen et al 



2009 



2009 



200 61 120081: iTreu et alj|2007t Uahnke et al. 
Greene et alj|2009t iMerloni et all I2009T iDecarli et al. 
Bennert et al.ll2009l) . With a few exceptions, most of these 
studies report a strong positive evolution in BH mass for fixed 
bulge properties, which reaches as high as ~ 0.6 dex offset in 
BH mass from the local relation at the highest redshift probed 
Cz~2). 

This observed strong evolution is difficult to understand in 
two aspects. For the M m - a relation, nume rical simulations 
based on self-regulated BH growth (e.g., Robertson et al 
2006) and other theoretical arguments (e.g., IShankar et al 
2009) favor a mild evolution in the normalization of the 
M. - a relation, in conflict with the otherwise claimed strong 
evolution. For the M, -Mb u i ge or M. -Lbuige relation, many 
of the observed hosts at earlier epoches are already bulge- 
dominated 2 , hence unless their bulges continue to grow a con- 
siderable amount over time, it is difficult to understand how 
these systems would have migrated to the local relation. 

An alternative explanation for this observed strong evolu- 
tion is that it is caused, at least partly, by the systematics in- 
volved in deriving both host properties and BH masses. All 
these studies mentioned above focus on broad line quasar 
samples, with host properties measured from the quasar-light 
subtracted images or spectra. The systematics in subtracting 

1 Hubble Fellow 

2 We note t hat there are also some observed hosts which are likely to be 
late-type (e.g., Sanchez et al 2004; Letawe et al 2007; Merloni et al. 2009). 



the quasar light, as well as in converting observables (such as 
photometric colors or spectral properties) to physical quanti- 
ties (such as bulge stellar mass) is a potential contamination 
to the final results. More importantly, the samples are pre- 
dominately selected in quasar luminosity and hence are not 
unbiased samples for studies of the BH scaling relations. 

One statistical bias resulting from using quasar-l uminosity 
selecte d samples was discussed in detail in lLauer et all 
d2007l) . Because there is cosmic scatter in the BH scaling 
relations and because the galaxy luminosity (or stellar mass, 
velocity dispersion, etc) function is bottom heavy, there are 
more smaller galaxies hosting the same mass BHs than the 
more massive galaxies. Hence when selecting samples based 
on quasar luminosities, low-mass BHs are under-represented, 
leading to an apparent excess in BH mass at fixed host prop- 
erties. If the cosmic scatter in BH scaling relations increases 
with redshift and reaches > 0.5 dex, it can in principle account 
for all the ex cess in BH mass observed for the mos t luminous 
quasars (e.g jLauer etaT]|2007l: IMerloni et al.ll2009h . 

A second statistical bias, resulting from the uncertainty 
in the BH mass estimates, was pointed out in Shen et al. 
(2008b, also see Kelly et al. 2009a). In all these stud- 
ies, the BH masses are estimated using the so-called virial 
method based on single-epoch spectra. In this method, 
one assumes that the broad line region (BLR) is virial- 
ized, and the BH mass can be estimated as M v ; r w G~ l RV 2 , 
where R is the BLR radius and V is the virial velocity; 
one further estimates the BLR radius using a correlation 
between R and the continuum luminosity L, i.e., R oc L Cl , 
found in local reverberation mapping (R M) samples (e.g., 
iKaspi et alJl2005l: iBentz etail 2006. 2009), and estimates the 
virial velocity using the width of the broad lines. In this 
way, one can estimate a virial BH mass using single-epoch 
spectra: log M vil = log R + 2 log(FWHM) + const = C\ log L + 
2 log(FWHM) + const. The coefficients are calibrated empiri- 
cally using RM samples or inter-calibrations between various 
lines (e.g.. iMcLure & Jarvisl l2002t IVestergaard & Peterson! 
l2006tlMcGilletaUl2008l) . " 
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Even though the virial method is widely used, these BH 
mass estimates based on several lines (usually Ha, H/3, 
Mgn and Civ) have a non-negligible uncertainty avir > 
0.4 dex, when compared to RM masses or masses derived 
from the M, — a r elation (e.g.. lVestergaard & Petersonll2006t 
lOnken et al.ll2004l ). This uncertainty must come from the im- 
perfectness of using luminosity and line width as proxies for 
the BLR radius and virial velocity, i.e., there are substan- 
tial uncorrelated variations ai in luminosity and variations 
ctfwhm m l me width, which together contribute to the overall 
uncertainty cr v i r , i.e., 



Cvir = V (C\ er L ) 2 + (2ct FW hm) 2 



(1) 



One can then imagine a statistical bias will arise if the un- 
derlying active black hole mass function (BHMF) is bottom 
heavy. In particular, the variations (uncompensated by the 
variations in FWHM) of luminosity at fixed true BH mass, 
ai, will scatter more lighter BHs into a luminosity bin than 
heavier BHs, and bias the mean BH mass in that bin. This 
is the Malmquist-ty pe bias (or Eddington bias) emphasized in 
IShen et aH(l2008bl sec 4.4), which has received little attention 
in the studies on the evolution of BH scaling relations to date. 

In this paper we examine the impact of this mass bias on the 
observed evolution in BH scaling relations with realistic mod- 
els for the underlying true BHMF and quasar luminosity func- 
ti on (LF). In j|2]we r eview the statistical mass bias discussed 
in lShen et al.1 (12008b) and demonstrate its effects with simple 
models; in §|3]we consider more realistic intrinsic BHMF and 
quasar LF, and estimate the mass bias as functions of lumi- 
nosity and redshift; we discuss the impact of this mass bias 
and conclude in §|4] Throughout the paper we use cosmologi- 
cal parameters f2o = 0.3, Da = 0.7 and Hq = 70 kms -1 Mpc -1 . 
Luminosities are in units of ergs" 1 and BH masses are in units 
of Mq. Unless otherwise specified, "luminosity" refers to the 
bolometric luminosity, and we are only concerned with the 
active BH population. 

2. THE BLACK HOLE MASS BIAS 

Denoting A = logL where L is quasar luminosity (bolomet- 
ric or in a specific band), m = logM. the true BH mass, and 
m e = logMyir the virial BH mass, the probability distribution 
of virial mass m e given true mass m is: 



p {) (m e \m) = (Ina^) 1/2 exp 



(m e -m) 2 



2a 1 

vir 



(2) 



where <7 v i r is the uncertainty of the virial BH mass. 

The probability distribution of true BH mass m given virial 
mass m e is: 



pi(m\m e ) =po{m e \m)^u{m)\ / p (m e \m)^ M (m)dm ^ 

(3) 

where 'I'm('w) = dn/dm is the true BHMF. If the sample is 
selected irrespective of luminosity (i.e., no flux limit) and so 
we can see all BHs with arbitrary m, then for a power-law 
true BHMF ^ M (m) oc 10 m7M and at fixed m e , the statistical 
bias between the expectation value (m) and virial mass m e is 
simply: 



AlogM. = m e -(m) = -ln(10)7 M er 2 



(4) 



In reality the sample is usually selected in quasar luminos- 
ity. If the distribution of luminosity at a given true BH mass 



for broad line quasars is g(\\m) where J g(X\m)dX = 1, the 
expectation value of m at fixed luminosity A is: 



(m)) 



J mg(X\m)^> M(m)dm 
J g(X\m)^J M(ni) dm ' 



(5) 



and the luminosity weighted average of m in the quasar sam- 
ple is: 



J^dX)dX 



(6) 



where Ai and A2 are the luminosity limits in the sample, and 
the quasar luminosity function ^ = dn/dX is 



^l(A)= / y M (m)g(X\m)dm. 



(7) 



On the other hand, the virial BH masses depend on lumi- 
nosities, and the sample averaged virial mass is: 



I>l(X) 



f m e f(m e I X)dm e 
j f(m e \X)dm e 



dX 



(8) 



where f(m e \X) is the probability distribution of m e given lu- 
minosity A and J ' f(m e \X) dm e = 1. In the absence of a theoret- 
ical model for f(m e \X), here we adopt an empirical recipe for 
f(m e I A): recall that the virial BH mass is expressed in terms of 
luminosity and FWHM. Neglecting the scatter from convert- 
ing bolometric luminosit y to continuum lumino sity, which is 
typically ^0.1 dex (e.g., Richards et al. 2006a), we have 



m e = C\X + 2 log(FWHM) + C 2 



(9) 



Here C\ « 0.5-0.7 is the slope in the measured luminosity- 
radius rel ation for BLRs and C2 is calibrated empiri- 
cally (e.g.. iKaspi et al.ll200l 120071: iBentz et alJl2006l l2009t 
iMcLure & Jarvisl 120021 : IVestergaard & Peterson! 120061 ). For 
statistical quasar samples, the distributio n of FWHMs does 
not depend on luminosity much (e.g., IShen et all [2008b; 
iFine et alJl2008T) . i.e., for any given luminosity, the FWHM 
distribution has a constant mean value and a log-normal scat- 
ter itfwhm around the mean. Therefore at fixed luminosity, 
the value of m e is independent of the true mass m, i.e., 



/(m e |A) = (27rai ne r 1/2 exp 



[m e -(dA + C 2 )] 2 
2a}- 



(10) 



where au ne = 2<7fwhm is the scatter in virial mass resulting 
from the scatter in FWHM (ctfwhm) at fixed luminosity, and 
the constant mean value of FWHM is absorbed in C2. 

Eqns. (O-© have analytical results for specific forms of 
^>M(m) and g(X\m). For a power-law true BH mass function 
\&Af(m) oc 10"' 7M , and a power-law relation between true mass 
and luminosity (A) = am+b with a constant log-normal scatter 
a L , i.e., 

g (A| OT ) = ( 27 ra 2 )-/ 2 exp(- [A - ( rr fo)]2 i, dD 



2ai 



Eqns. ©-© yield 



* L (A) oc 10 A7M/fl ; 



(12) 



(m) = 



7Mln(10)er 2 -a& 



1 



7M ln(10) 

A 2 io A ^*/ fl -AilO Al7M / a 
a(10 A27M / a - 10 Al7 "/°) 



;(13) 
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(m e ) =C 2 - 



C\CL 



+ 



7Mln(10) 

Ci(A 2 10 A ^"/ fl -AilO Al T*/ fl ) 



(14) 



Note throughout this paper we assume this constant scatter in 
luminosity, ox, is solely the variations that are uncompensated 
by the corresponding variations in FWHM, i.e., it is one of the 
two sources of the virial uncertainty er v j r (Eqn.[TP. 

There are eight parameters that determine the bias 
AlogM, = (m e ) - (m): ox, jm, a and b are for the underlying 
true BH mass function and luminosity (or Eddington ratio) 
distributions at fixed true mass, which are model assumptions 
and must be tuned to produce the observed quasar luminosity 
function; C\ and C2 are determined from the empirically cali- 
brated virial estimators and measured FWHM distributions in 
statistical quasar samples; Ai and A2 are observational win- 
dows determined from the specific quasar sample under in- 
vestigation. A further constraint comes from the assumption 
that virial mass is an unbiased estimator of the true BH mass 
(e.g., Eqn.[2]i, which requires (e.g., consulting Eqns.|2l[T0land 



C\ = \/a, C2 = —b/a 



Thus the bias reduces to 



AlogM. = -7 M ln(10)cr£/a 2 



(15) 



(16) 



from E qns. ( fT3] l and (fl4l (see also sec 4.4 in IShen et aT] 
2008b), which is independent on luminosity. 

For demonstration purposes here we take a model with 
parameters ctl = 0.4, 7^ = -4, a = 2 and b = 28.4. These 
parameters produce a slope of 72, = -2 in t he quasar LF, 
close to the observed bri ght-end slope (e.g., iRichards et al.l 
2006b; Ho pkins et al.l2007l) . and constants C\ = \ja= 0.5 and 
C2 = -b/a = -14.2, which are consistent with the commonly 
used virial m ass calibrations and the observed FWHM distri- 
butions (e.gjMcLure & Jarvisll2002t Vestergaard & Peterson! 
120061 iMcGill et alj|2008t IShen et al.ll2008bl) . With these pa- 
rameters, the difference between the sample-averaged true 
and virial BH masses is (m e ) - (m) ra 0.37 dex. 

The derivation of this virial mass bias is independent from 
the derivation of the bias discussed in lLauer et al.l (120071) . To 
illustrate this, let us consider again simple power-law dis- 
tributions and Gaussian (lognormal) scatter. We start from 
a power-law distribution of quasar host galaxy property s 
(s = logMbuige for instance), 



The distribution of true BH mass m at fixed s is: 

[m-(C 3 5 + C 4 )] 2 



/7 (m|5) = (2w 2 r 1/2 exp 



(17) 



(18) 



where the BH scaling relation is (m) = C3S+C4 with a cosmic 
scatter a^, where C3 and C4 are constants. Given the power- 
law distribution of s (Eqn. [17), the distribution of true BH 
mass m is then also a power-law v I'm('«) oc 10™ 7m where 7m = 
75/C3. Assuming the distribution of luminosity A at fixed m is 

3 It is possible to incorporate an additional correlated variation term in 
luminosity at fixed true mass, which is compensated by the variations in 
FWHM and hence does not add to the virial uncertainty, but this will makes 
the mass bias discussed in 0even more severe (see later discussions). 



given by Eqn. ( fTTT i we can derive the distribution of A at fixed 

s as: 

A>(A|s)= / g(X\m)p (m\s)dm 



■■[2TT(al + a 2 a 2 )] 1/2 exp<' 



2(a 2 + fl 2 «T 2 ) 



Therefore the distribution of s at fixed luminosity A is 



(19) 



/>i(i|A) = po(A|i)* s (i) / po(X\s)^ s (s)ds 



2\-i/2 I [i-(C 5 A+C 6 )] 2 



= (271-ofr ' exp 



2a 2 



, (20) 



where 



a 2 + a 2 a 2 , 1 



a i = L iri ^ C 5 = —,C 6 = ln(10) 7s crf 



2 b + aCa, 
aCi 



(21) 



Similarly the distribution of m at fixed A is: 

pi(m\X) = g(X\m)^ M (m)( / g(X\m)^ M (m)dm 



where 



= (27rcr|r 1/2 exp 



— k , Oj , C8 : 

a a 



[m-(C 7 A + C 8 )] ; 
2a\ 



ln(10) 7s cr 2 b 



a 2 Ci 



, (22) 



(23) 



Hence the Lauer et al. bias, i.e., the excess in true BH mass at 
fixed luminosity is simply: 

A logM., Lauer = (m) I A - (C 3 (s) I a + CO 

= -ln(10) 7s cT 2 /C 3 . (24) 

Since 75 is negative, the Lauer et al. bias states that at fixed 
luminosity, the average true BH mass is already biased high 
with respect to the expectation from the mean BH scaling re- 
lation. On top of that, the virial mass bias (Eqn. [Tol l states 
that the average BH mass estimate is further biased high with 
respect to the true mass. 

3. MORE REALISTIC MODELS 

The simple estimate of the BH mass bias discussed above 
neglects the curvature in the intrinsic BHMF and quasar LF, 
hence is only valid at the bright end of the LF. In particu- 
lar since the quasar LF flattens below the break luminosity 
Abreak ~ 46, we expect that the mass bias becomes less se- 
vere towards fainter luminosities. It is important to realize 
that in essentially all of the studies mentioned in ^T] larger 
offsets usually occur at higher redshift, when their samples 
are sampling the higher-mass end of the intrinsic BHMF. It is 
conceivable then that a false evolution will arise simply from 
the increasing mass bias towards the high-mass tail. 

To estimate the BH mass bias at various redshifts and lumi- 
nosity sampling ranges, we assume an underlying BHMF and 
a model for the luminosity distribution at fixed true BH mass 4 , 

4 We note that in principle one should use an underlying BHMF that 
is constrained by the measured BHMF based on virial masses (see, e.g., 
Kelly et al. 2009a b). Unfortunately, the derived BHMF is poorly constrained 
at M, < 10 s Mq for current flux-limited quasar samples. 
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TABLE 1 
Model Parameters 
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several redshifts (green lines). The points are the d erived bolomet ric LF data 
based on optical and soft X-ray data compiled by Ho pkins et a n j2007l and 
references therein). 
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FIG. 2. — Estimated BH mass bias using our best-fit models of the true 
BHMF and luminosity distributions, for luminosity-limited quasar samples 
with A > Ao = logLo at various redshifts. Shown here are the results for 
three fixed values of ctl = 0.3, 0.4 and 0.5 dex. The bias increases when 
the luminosity or 07, increases. We also plo t the observed BH mass ex- 
cess at fixed host properties from Woo et al. (2006, 2008, filled square), 
Salviander et al. (2007, open circles), Merloni et al. 1 2009, filled trianges), 
and Decarli et al. (2009, open squares) at the corresponding redshifts. Error 
bars are the standa rd deviations of objects in their samples. Note that the 
data from Salviander et al. (2007) are already binned, hence the dispersion is 
substantially smaller. 

using the observed quasar LF as constraints. This forward- 
modeling procedu re is similar to the sim ulations done in 
IShen etal.l (l2008b1) and lKellv et all J2009allbh . 

We start with a broken power-law model for the true 
BHMF: 



* M (m) 



DC 



(25) 

\Q-lM\(m-m„) + JQ-7m2(«i-'«*) ' 

where 7^1 and 7^2 are the low-mass and high-mass end slope, 
and m* is the break BH mass. For the luminosity distribution 
g(X\m) we still use the Gaussian form in Eqn. ( fTTT i. For sim- 
plicity, we assume that the set of six parameters, 7^1, Jmi, 
to*, a, b and <tl, are only functions of redshift. At a fixed 
redshift, a set of the six parameters determines the shape of 
the quasar LF (0, which is to be constrained by observa- 
tions. Assuming that the commonly used virial BH estima- 
tors are unbiased, we fix a = 2 and b = 28.4, which corre- 
sponds to a BLR radius-luminosity relation R oc L° 5 and virial 
coefficients and FWHM distributions consistent with previ- 
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NOTE. — The x 2 fits are done at individual re dshifts against the bolo- 
metric LF data compiled in Hopkins et al. ( 2007), i.e., we do not fit a 
global broken power-law BHMF model. For certain values of ul and 
redshifts, the high-mass end slope is poorly constrained, in which case 
we fixed the high-mass end slope to be 7^2 = —6. 



ous work (e.g.. IShen et ai]|2008bl: iFine et al.ll2008l) . We also 
fix the value of 07, during model fitting. Therefore there are 
four free parameters: 7^1, 7m2, m*, and the normalization of 
the BHMF, $ (in units of Mpc" 3 dex _1 ). We list the best-fit 
model parameters in Table 1 for different values of <jl and at 
several redshifts. We note that these models are mainly used 
to demonstrate the effects of the virial mass bias, and should 
not be interpreted as our attempt to constrain the true active 
BHMF. 

Fig.Q]shows the model LF with fixed <tj, = 0. 4 and compar- 
ison w ith the bolometric LF data compiled in Hopkins et al. 
d2007l and references therein), at several redshifts. At each 
redshift, we determine 7^1, 7m2, tn* and $0 by minimizing 
the x 2 between the model and the data. As discussed earlier, 
this scatter tj£ refers to the variations that are uncompensated 
by the corresponding variations in FWHM. If we were to in- 
corporate an additional correlated variation term in luminos- 
ity at fixed true mass, which are compensated by an additional 
variation in FWHM, we would need a steeper high-mass end 
slope for the true BHMF in order not to overshoot the bright- 
end LF, which then makes the mass bias even worse. Also, 
the amount of this additional broadening i s limited by the al- 
ready narrow dis tributions of FWHM (e.g. JShen et al.f2 008b; 
IFine et ai1l2008l) and the requirement that cr v i r > 0.4. 

Given the model true BHMF and luminosity distribution, 
we can now estimate the bias A logM. = (m e )- (to) as a func- 
tion of the sampled luminosity range using Eqns. © and 
(0. Fig. [2] shows AlogM. for luminosity-limited quasar 
samples with A > Ao = logLo at several redshifts and for 
fixed <jl = 0.3, 0.4, 0.5 dex. The bias generally rises when 
luminosity increases and approaches the asymptotic value 
- 7M2ln(10)<7^/a 2 at the bright end. Because we are fitting our 
model against the bolometric LF data at individual redshifts, 
we do not expect identical results since there will be system- 
atics involved in deriving the bolometric LF data from band 
conversions and from the assumed quasar spe ctral energy dis- 
tribution (see details in Hopkins et al. 2007). Nevertheless, 
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this analysis demonstrates that the mass bias (m e ) - (m) could 
be substantial for bright luminosities and for large dispersion 
<?£ in the luminosity distribution at fixed true mass. 

4. DISCUSSION 

The virial mass bias discussed here depends on the ampli- 
tude of the scatter c>l. This is the variation in luminosity at 
fixed true BH mass which is not compensatedby the variation 
in FWHM. How large is a L l For the best studied local RM 
samp l es and for H/3 only (e.g.. lKaspi et al.ll2005HBentz et al.1 
120061 l2009t) . the scatter of luminosity around fixed BLR size 
is on the order of < 0.2 dex. However, this is for the best 
cases. We generally expect larger ctl for statistical quasar 
samples with single-epoch spectra for the following reasons: 
1) single-epoch spectra have larger scatter than averaged spec- 
tra in RM samples; 2) these samples usually cover a luminos- 
ity range that is poorly sampled by current RM data; 3) the 
UV regime of the spectrum (in particular for Civ) has more 
peculiarities and larger scatter in the R-L relation than the H/3 
regime, but is usually used at high redshift due to the drop off 
of H/3 in the optical spectrum; 4) even if some variations in lu- 
minosity truthfully trace the variations in BLR radius R, they 
may still be uncompensated by FWHM if, for instance, there 
is a non-virial component in the broad line which does not 
respond to BLR radius variatio ns (especially for Civ, see dis- 
cussions in IShen et aLll20 08b): 5) in many statistical quasar 
samples, the spectra have low signal-to-noise ratios, which 
lead to increased uncorrelated scatter between the measured 
luminosity and FWHM, and bias the virial mass estimates 
more. Hence it is very plausible that <7l > 0.4 in most cases, 
which then accounts for > 0.2 dex in cr v i r (e.g., Eqn.[TJ. 

Therefore we conclude that the mass bias discussed here 



contributes at least 0.2-0.3 dex at Lboi ^ 10 46 ergs" 1 , com- 
parable to the statistical bias resultin g from the cosmic scat- 
ter in the BH scaling relations (e.g., lLauer e t al. 2007]). In 
Fig. |2] we also show the "observed" BH mass offset (excess 
at fixed galaxy properties) from the literature. We note that 
the amount of the uncorrelated scatter oi might increase with 
redshift both due to the switch from the H/3 line to the more 
problematic UV lines and due to the often decreased spectral 
quality at high redshift. Since the virial mass bias is inde- 
pendent of the Lauer et al. bias, the two statistical biases to- 
gether can account for a large (or even full) portion of the BH 
mass excess seen in the data, hence no exotic scenarios are 
needed to explain this strong redshift evolution. It is possible, 
however, that there is still a mild evolution in these BH scal- 
ing relations, which is inherent to the cosmic co-formation 
of SMBHs and their hosts. But given these statistical biases, 
and given other systematics with singl e-epoch virial estim a- 
tors (not discussed here, but see, e.g.. lDennev et a l. 2009), it 
is premature to claim a strong positive evolution in the BH 
scaling relations. 
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